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Abstract 

We prove that the spectral gap of the Swendsen-Wang process for the 
Potts model on graphs with bounded degree is bounded from below by some 
constant times the spectral gap of any single-spin dynamics. This implies 
rapid mixing for the two-dimensional Potts model at all temperatures above 
the critical one, as well as rapid mixing at the critical temperature for the 
Ising model. After this we introduce a modified version of the Swendsen-Wang 
algorithm for planar graphs and prove rapid mixing for the two-dimensional 
Potts models at all non-critical temperatures. 



1 Introduction 

We study the mixing properties (in terms of the spectral gap) of the Swendsen- 
Wang dynamics for the g-state Potts model on graphs with bounded degree. The 
Swendsen-Wang process is probably the most widely used algorithm to generate 
(approximate) samples from the Potts model because, empirically, it seems to be 
rapid for all graphs and temperatures, at least for small q. But this is in general not 
true. It is known that the process does not mix rapidly at some critical temperature 
on the complete graph for all q > 3, see Gore and Jerrum [GJ97] . and on rectangular 
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subsets of the hypercubic lattice Z d , d > 2, for q sufficiently large, see the recent 
paper of Borgs, Chayes and Tetali |BCT10j . Furthermore, Li and Sokal (LS89J 
proved that there is a lower bound on the mixing time in terms of the specific 
heat. Results on rapid mixing of the Swendsen-Wang dynamics are also rare. Rapid 
mixing is proven at all temperatures on trees and cycles for any q, see Cooper and 
Frieze |CF98j . on narrow grids for any q, see Cooper et al. |CDFR00~] . and on the 
complete graph for q = 2, sec [CDFROOJ and Long et al. [LNP07J. For rapid mixing 
results for sufficiently high or low temperatures see e.g. Martinelli |Mar92j . Cooper 
and Frieze |CF98j and Huber |Hub03j . 

In this article we prove that the spectral gap of the Swendsen-Wang dynamics 
on graphs with bounded degree is bounded from below by some constant times the 
spectral gap of the (single-site) heat-bath dynamics. The main result of this paper 
can be stated as follows, see Theorem [TJ 

Suppose that P (resp. Phb) is the transition matrix of the Swendsen-Wang (resp. 
heat-bath) dynamics, which is reversible with respect to the Boltzmann distribution 
7Tp q of the g-state Potts model on the graph G at inverse temperature (3. Then there 
exists a constant c > 0, depending only on q and (3, such that 

A(P) > c A A(P HB ), 

where A(-) denotes the spectral gap and A is the maximal degree of the underlying 
graph G. In Theorem [TJ (Subsection 4.3.2) we prove this statement also for the case 
that we have only a single vertex v* of large degree, where A is replaced by the 
second largest degree and Phb is replaced by the transition matrix of the heat bath 
chain that leaves the color at v* constant. 

In particular, this implies (see Corollary [2]) that the Swendsen-Wang process is 
rapidly mixing above the critical temperature for the two-dimensional Potts model 
(q > 2) and at the critical temperature for q = 2. Furthermore we get rapid mixing 
for graphs with maximum degree A if p — 1 — e^ 13 < 1/A. If the graph is planar, 
then p = 1 - e~P < 1/(3^(^-3)) is sufficient for A > 6 (see Corollary ED . This 
improves over the 1/(3 (A — 1)) bound of [Hub03]. 

One may believe (and numerical results suggest) that there could be a sharper 
bound on the spectral gap of the Swendsen-Wang with an additional factor N (= 
number of vertices of the graph) on the right hand side, i.e. A(P) = ^(A^A(Phb))- 
But proving this does not seem to be possible with the techniques of this paper. 

We conjecture that the spectral gap of Swendsen-Wang is (almost) the same 
at high and low temperatures, but we were not able to prove this. Nevertheless 
we present a modified version of the Swendsen-Wang algorithm for planar graphs 
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that makes an additional step on the dual lattice. We prove that this process is 
rapidly mixing for the two-dimensional Ising model at all temperatures and for the 
two-dimensional Potts model for all non-critical temperatures (see Corollary [TTj) . 



2 Preleminaries 

In this section we introduce the models and notations to analyze the Swendsen- 
Wang algorithm. First we define the Potts (Ising) and the random-cluster model 
and explain their connection. 

2.1 The models 

Fix a natural number q > 2 and a graph G = (V, E) with finite vertex set V and 
edge set E C ( 2 ), where ( 2 ) is the set of all subsets of V with 2 elements. Let 
N := |V|. As an example one can have in mind the two-dimensional square lattice 
G L of side length L, i.e. V = {1, . . . , L} 2 C Z 2 and E = {{u, v} G V 2 : \u - v\ = 1}, 
where | • | denotes the l\ norm. 

If we consider more than one graph, we denote by Vq (resp. Eq) the vertex (resp. 
edge) set of the graph G. Furthermore we write u -f* v if u and v are neighbors in 
G, i.e. {u,v} G E, and by deg G (t>) we denote the degree of the vertex v in G. The 
q-state Potts model on G is defined as the set of possible configurations flp = [q] v , 
where [q] := {1, . . . , q} is the set of colors (or spins), together with the probability 
measure 

7^0) := ^ q (a) = ^ ^ exp <j ^ 1 (a(u) = a(v)) 



u,v: u+->V 



for a G flp, where Z is the normalization constant (also called partition function) 
and (3 > is called the inverse temperature. This measure is called the Boltzmann 
(or Gibbs) distribution. If q = 2 the Potts model is called the Ising model. 

A closely related model is the random cluster model (also known as the FK- 
model), see Fortuin and Kasteleyn |FK72j . It is defined on the graph G = (V,E) 
by its state space fl^c — {A : A C E} and the RC measure 



\A\ 
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where p G (0, 1), C(A) is the number of connected components in the graph (V, A), 
counting isolated vertices as a component, and Z(-,-,-) is the same normalization 
constant as for the Potts model (see [Gri06j Th. 1.10]). For a detailed introduction 
and related topics see |Gri06] . 

There is a tight connection between the Potts model and the random cluster 
model. Namely, if we set p = 1 — e - ^, we can translate a Potts configuration a ~ irp 
to a random cluster state A ~ fi p and vice versa. To get a Potts configuration a G flp 
from A G £7rc assign a random color independently to each connected component 
of (V, A). For the reverse way include all edges e = {ei, t<{\ G E with cr(ei) = cr(e2) 
to A with probability p. Hence sampling a Potts configuration according to 7Tg is 
equivalent to sampling a RC state from \x v if both models are defined on the same 
graph G and p = 1 — e _/3 . For a proof see [E588j. 



2.2 The heat-bath dynamics 

The heat bath dynamics is probably the most studied Markov chain related to the 
Potts model (especially for q = 2). This is because its mixing time is related to some 
properties of the underlying model. See e.g. the monograph of Martinelli |Mar99j . 
The transition matrix of the heat-bath chain on Qp is defined by 

P UB (a,a ) .= P nBAq (a,a ) = - (1) 

where a v,k {v) — k G [q] and a v,k {u) = a(u), u^v. Otherwise Phb(o", t) = 0. 

This transition matrix describes the process that, at each step, chooses one vertex 
of the graph G uniformly at random and changes only the color of this vertex with 
respect to the conditional distribution of 7Tg, given the color of all other vertices. It 
is easy to prove that this Markov chain is reversible with respect to n^. 

The spectral gap of Phb is well known on some classes of graphs. For positive 
results see e.g. |Mar99j . |BKMP03j . |MT10j and jLSlO] . We will state a few results 
in Section [3l 



2.3 The Swendsen-Wang algorithm 

Now we turn to the Swendsen-Wang algorithm [SW87]. First we state the coupling 
of the Boltzmann distribution tt^ q and the random-cluster measure of Edwards 
and Sokal [ES88J. Let us define 

Q(A) := [a G fip : a(u) = a(v) V{u,v} G A}, A G fi RC , 
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and 

E ( a ) '■= {{u,v} <E E : a(u) = a(v)}, cr G fi P . 

Obviously, we have for a G Qp and A C E that cr<EQ(A) <^ ic£(c). Let a G H P , 
A G f^Rc an d p = I — e~P, then the joint measure of (a, i) 6 fl P x f2 RC is 

p 



The marginal distributions of p are exactly Tip and /x p , respectively. The Swendsen- 
Wang algorithm is based on this connection of the random cluster and Potts models 
and performs the following two steps. 

1) Given a Potts configuration o~ t G Qp on G, delete each edge of E(o~ t ) indepen- 
dently with probability 1 — p — . This gives A G JIrc- 

2) Assign a random color independently to each connected component of (V, A). 
Vertices of the same component get the same color. This gives cr t +i £ ^p- 

This can be seen as first choosing A with respect to the conditional probability 
of p given a t and then choosing a t +i with respect to p given A. If we define the 
|Qp| x |f2Rch m atrix 

Ta^A) := = p M (1-^)1^)1-1^1 1{A c (2 ) 

and the |J1rc| x I^p (-matrix 

T ^ (Aa) := = <i~ c{A) U°en{A)), (3) 

then the transition matrix of the Swendsen-Wang dynamics (on VLp) is defined by 
P(a,r) := P^r) = T G;P;q T^ q (a,r), a,ren P . (4) 
Note that the transition matrix of the Swendsen-Wang dynamics on Qrc is given by 
P(A,B) := ?,(A5) = T^ q T G , p , q (A,B), A,Ben RC . (5) 
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2.4 Spectral gap 



In the following we want to estimate the spectral gap of certain transition matrices of 
Markov chains. For an introduction to Markov chains and techniques to bound the 
convergence rate to the stationary distribution, see e.g. [LPW09J. If the transition 
matrix P with state space Q is ergodic, i.e. irreducible and aperiodic, and reversible 
with respect to tt, i.e. 

n(x) P(x, y) = ir(y)P(y,x) for all x, y E Q, 

we know that 1 = £o > £1 > • • • > £|fi|-i > — 1, where the & are the (real) eigenvalues 
of P. The spectral gap of the Markov chain is defined by X(P) = 1— max{£i, K|n|-i }• 
If we are considering simultaneously a family of graphs G = (Vq, Eg), we say that 
the chain is rapidly mixing for the given family if A(P) -1 = 0(|Vg| ) for some 
C < oo. 

The eigenvalues of the Markov chain can be expressed in terms of norms of the 
operator P that maps from L 2 (7r) := (R n , \\ ■ \\ n ) to L 2 (vr), where scalar product 
and norm are given by = J2 x£n f(x)g(x)ir(x) and \\ f\\l := J2 xen f(x) 2 7r(x), 

respectively. The operator is defined by 

pf{x) ■.= J2 p ( x iy)f(y) 

yen 

and represents the expected value of the function / after one step of the Markov 
chain starting in x 6 fi. The operator norm of P is 

\\P\\ n := \\P\\l 2 (^l 2 (7t) = max \\Pf\U 

and we use || • ||„. interchangeably for functions and operators, because it will be clear 
from the context which norm is used. It is well known that X(P) = 1 — \\P — S^W 
for reversible P, where S n (x,y) = 7t(y). We know that reversible P are self-adjoint, 
i.e. P = P*, where P* is the (adjoint) operator that satisfies (f,Pg) n = (P*f,g) n 
for all /, g E L 2 (it). Note that if we look at Ta lPl q from (j2J) as an operator mapping 
from L 2 (fi p ) to L 2 (vr^), then T G is the adjoint operator of T Gp (? . This proves that 

P and P from and (jSJ) have the same spectral gap. 

3 Main result 

In this section we prove by comparison that the spectral gap of the Swendsen-Wang 
dynamics is bounded from below by some constant times the spectral gap of the 
heat-bath chain. Therefore we fix some value of f3 > and q > 2. 
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We will prove the following theorem. 

Theorem 1. Suppose that P (resp. Pub) is the transition matrix of the Swendsen- 
Wang (resp. heat-bath) dynamics, which is reversible with respect to k9 . Then 

A(P) > c sw A(Phb), 

where 

csw = csw(G,p,q) := ^ (qe w ) 

with 

A := max deg G (v). 

veVc 

By this result we get that the Swendsen-Wang dynamics is rapidly mixing (on 
graphs with bounded degree), if the usual heat-bath dynamics is rapidly mixing. 
Unfortunately, this result seems to be off by a factor of N, because we compare the 
SW dynamics with a Markov chain that changes only the color of one vertex of the 
graph per step. A similar comparison with the systematic scan heat-bath dynamics 
would give a better bound (see e.g. |Hay06| ) . But this does not seem to be possible 
with our techniques. 

Now we will state some results that follow from Theorem [U 

The first corollary deals with the two-dimensional Potts model. Recall that the 
two-dimensional square lattice Gl of side length L is given by V = {1, ... , L} 2 C Z 2 
and E = {{u, v } e V 2 : \u — v\ — 1}. 

Using the known lower bounds on the spectral gap of the heat bath dynamics 
on the square lattice [MOS941 IBD10[ ILSlOj . we can obtain the following result: 

Corollary 2. Let Gl be the square lattice of side-length L, N = L 2 . Then there 
exist constants Cp = cp(/3, q) > and m > such that 

• KPp, L q T l < c?N for (3 < c (q) = ln(l + Jq) 

• K p p, L q )~ l < cN m forq = 2 and (3 = (3 C (2), 

where c" 1 = c S w(Gl, c , 2). 

Not much is known about the constant eg, but there is an explicit bound on the 
exponent m in terms of an (unknown) crossing probability in the random cluster 
model, see |LS10j . 
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Proof of Corollary [H We have to show that the heat-bath dynamics is rapidly mix- 
ing in the cases stated above. The result for q = 2 at the critical temperature is 
given by Lubetzky and Sly [LSlOj . For the high temperature result we need the 
recent result of Beffara and Duminil-Copin [BDTUJ Th. 2]. They prove that one 
has exponential decay of connectivity in the RC model on Z 2 for all q > 1 and 
p < = i _ e -Mi), Together with |Ale98l Th. 3.6] and |MUS94l Th. 3.2] we get 

that there exists a constant cp > such that A(Phb)" 1 < c@ N on Gl for /3 < (3 c (q). 
Thus the result follows from Theorem [U □ 

The next corollary relies on a result of Hayes Hay06|, who gives a simple condi- 
tion on (3 for rapid mixing of the heat-bath dynamics for the Ising model on graphs 
of bounded degree. Especially, we use his result for planar graphs. 

Corollary 3. Let G be a graph with maximum degree A. Then for every q > 2 and 
e > the spectral gap of the Swendsen- Wang dynamics satisfies 

->G \-l - c 



A(P^)- 1 < -NlogN, 



where c 1 = \ c S w(G, /?, q), if 
or if G is planar, A > 6 and 



< 



1-e 



X 3(A - 3) 

If we state the bounds on the temperature in terms of the RC parameter p 



1— e P, this leads to the bounds p < (1 — e)/ A in general and p < (1 — e)/(3vA — 3) 
for planar graphs. 



Proof of Corollary^ In |Hay06 Prop. 14, Coro. 19] the result is given for q = 2 



and f3 < (1 — s)/p in terms of mixing time, where p is the principal eigenvalue of 
the graph. The bounds on p are stated therein. Note that the different definition 
of Tip causes the additional factor 2. To generalize this result to the g-state Potts 
model one only has to prove an inequality like in Observation 11 in |Hay06|. But 
this is easily done, e.g. by induction over q. Finally we have to bound the spectral 
gap in terms of the mixing time. For this see e.g. [LPW091 Th. 12.4]. □ 
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Remark 4. Note that one can read the statement of Theorem [T] also as an upper 
bound on the spectral gap of the heat bath dynamics. See e.g. |BCT10j for a result 
on slow mixing of the Swendsen-Wang dynamics. 



3.1 Proof of Theorem [T] 

See Diaconis and Saloff-Costc [DSC93J for an introduction to comparison techniques 
for Markov chains. For the comparison with the Swendsen-Wang dynamics we will 
analyze the Markov chain with transition matrix 

Q = P UB PP nB . (6) 

Since P BB and P are reversible with respect to 7Tg, Q is also reversible. 



Lemma 5. With the definitions from above we get 

HQ) > MPhb) > A(Phb). 

To prove a lower bound on the spectral gap of P it remains to prove X(P) > cX{Q) 
for some c > 0. For this we need an estimate of the transition probabilities of the 
Swendsen-Wang dynamics on G with respect to some subgraph of G. Therefore we 
prove the following lemma. 



G with E C E. Define P G := P9 Then 



Lemma 6. Let G = (V, E) be a graph and Gq = (V, E ) be a spanning subgraph of 

^G 

P,r 

c^ oI P g >,t) < P G (a,r) < ^ Eol P Go (a,r) 



for all a, t e Qp, where 

ci = ci(/?) := e- p 

and 

c 2 = c 2 ((3,q) := l + g ( e ^-l). 

Proof. The first inequality is already known from the proof of Lemma 3.3 in |BCT10j , 
but we state it here for completeness. Let p — 1 — e _/3 , then 

P G (a,r) = P [A[ (1 ~ p) m ° )HA[ q-° {A} t(A C E(a) n E(t)) 
Ace 

> P lAl ( x - p) lE((T)HAl q~ c(A) t(Ac E{a) n E{t)) 

AcE 

> (l- p )W«)\-\z°(°)\p Go (a,T) > (1-p)\ e \ e °\P Go (ct,t). 
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For the second inequality suppose for now E = E \ {e} for some e G E and note 
that C(A U {e}) > C{A) - 1. We get 

P G (a,r) = J2 p |A| (l-p) |£MHA| g- C{A) 

AcE(a)nE(r) 



AaE(o)C\E(r): 



\\E(*)\-\A\ -C(A) 



AcE(a)nE(r): 
e<£A 



< p |A'U{e}| ^ _ p yE(a)\-\A'u{e}\ q -C(A'u{e}) 

A'dE {a)f\E (T) 

+ (1 - p )\ E W\-\ A '\ q -°( A ') 



A'cEo(a-)nE (T) 



< 



qp ^ p lA ' 1 (1 - p )\ E ^\-\ A '\ q-°^ 

1 P A'CE (a)nE (r) 

+ (i-p) |BWH£ ° (ct)| p G oM 

(y^ + l) Pgo(^t) = (1 + q(e^-l))P Go (a,r). 



For |£'\£'o| > 1 one can iterate this technique \E\E \ times. This proves the 
statement. 

□ 

We use this lemma to prove that the transition probability from a to r is similar 
to the probability of going from a neighbor of a to a neighbor of r. Recall that cr v,k 
is defined by a v ' k (v) — k G [q] and a v ' k (u) = cr(u), u ^ v. 

Lemma 7. Let a,r G Hp, v G V G and k, I G [q]. Then 

Pg{<J v '\t v > 1 ) < cf SG{v) P G (a,r) 

with 

c 3 = c 3 ((3,q) := qe^-(q-l)^. 
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Proof. Define E v := {e e E G : v e e} and G v := (V,E \ E v ). Then, v G Vq is 
an isolated vertex in G„. By the definition of the Swendsen-Wang dynamics we get 
that 

Pg v (v v '\t v > 1 ) = PgMt). 
If we set E = E \ E v we get \E\E \ = deg G (t> ) and by Lemma |6] 

P G (a v > k ,r v > 1 ) < cf 8Giv) P Gv (o- V > k ,T"> 1 ) = cf SGiv) P G >,r) 
with Ci and Ci from Lemma [6j 

□ 

Now we are able to prove the main result. 

Proof of TheoremUl Because of Lemma [5]we only have to prove A(P) > csw A(Q). 
Let 

c := max — -, 

o-i,CT2,Ti,T 2 eOp -r(CT2,T2) 

where a ~ r X^ev l <T (' u ) ~~ r ( t ')l — 1- Note that Phb(c, t) 7^ if and only if 
cr ~ r. We get for a±, T\ G Qp that 

<5((Ti,ri) = 2j -PHB(0-l,Cr2)-P(0"2,r 2 )PHB(T2,Ti) 

f2,T2enp 

< cP((Ti,Ti) PhbOi,<X 2 ) Phb(t2,Ti) 

(7 2 ~(Tl T2~Tl 

< gcP(cri,ri). 

It is well-known that 

A(p2) = min {tr77r Var(/) ^° 

where Var(/) = ||/ - 3 Vp f\\l p and 

<7,T 
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We use P 2 instead of P, because this representation holds only if A(P) = 1 — £i(P) 
|DSC96j . From the calculation above we have £pz(f) > -^s £q 2 (/) f° r every / G 
L 2 (7r / 3) and so A(P 2 ) > -^X(Q 2 ). Since we have for reversible P that A(P) < 
A(P 2 ) < 2A(P), we get 

1 

2^c 2 

It remains to bound c. With C3 from Lemma [7] we get for cxi, cr 2 , 7~i, r 2 G fip with 
o"i ~ o"2 and n ~ r 2 



> ^a(q). 



P(<ruTi) ^ a-P^jTi) . 2A P(a 2 ,r 2 ) 



< ct „ < c: 



,2A 



P(^,r 2 ) ~ 3 P(<t 2 ,t 2 ) ^ P(a 2 ,r 2 ) 
Finally, 

w > > ^ c 3 4AA ^) ^ ^(^ 2/3 r 4A A(Q). 

This completes the proof. 

□ 



4 The Modified Swendsen-Wang 

With the techniques from the last section we do not get any positive result for low 
temperatures. However, the Swendsen-Wang dynamics seems to be rapidly mixing 
also at low enough temperatures. Unfortunately, we were not able to prove it. 

In this section we introduce a modified version of the Swendsen-Wang dynamics 
for planar graphs that is rapidly mixing for the two-dimensional Ising model at all 
temperatures. Roughly speaking the chain makes one step at high and one step at 
low temperatures. 

4.1 Dual graphs 

A graph G is called planar if one can embed it into IR 2 such that two edges of G 
intersect only at a common endvertex. We fix such an embedding for G. Then we 
define the dual graph Go = (Vd,Ed) of G as follows. Place a dual vertex in each 
face (including the infinite outer one) of the graph G and connect 2 vertices by the 
dual edge if and only if the corresponding faces of G share the boundary edge e 
(see e.g. |GrilOl Section 8.5]). It is clear, that the number of vertices can differ in 
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the dual graph, but we have the same number of edges. 

Additionally we define a dual RC configuration C Ejj in to a RC state 
A C £ in G by 

e G A <<=^ e D £ A D , 

where cd is the edge in Ed that intersects e in our (fixed) embedding. (By con- 
struction, this edge is unique.) See Figured] for the graph Gl with L = 3 and its 
dual graph (Gl)d together with 2 corresponding RC states. 
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Figure 1: Left: The graph G3 (solid) and its dual (dashed). 

Right: A RC state on G 3 (solid) and its dual configuration (dashed). 



It is easy to obtain (see |Gri06t p. 134]) that the random cluster models on the 
(finite) graphs G and Gr> are related by the equality 

^ q (A) = li%{A d ), (7) 

where the dual parameter p* satisfies 

1 — p* p 

The self-dual point of this relation is given by p s d(q) = i+^j which corresponds by 
p = 1 — e _/3 to the critical temperature of the g-state Potts model (3 c (q) = log(H-y^) 
on Z 2 |BD10j . In the following we denote by 0* the "dual" value of /3, i.e. we have 
p* = 1 — e _/3 *, p — 1 — e _/3 and (jHJ). If G is fixed, we write jx* instead of ll Gd . 
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4.2 The algorithm 

To describe the transition matrix of the modified Swendsen-Wang algorithm, we 
define the (square) 2' E ' x 2' £D '-matrix 

D(A,B) := 1(B = A D ), AcE,BcE D . 

By (j7j) we get that fi p ■ D = fi*,, respectively ■ D* = fi p (distributions are seen 
as row vectors), where D*(B, A) = D(A, B), A C E, B C E D . 

Now we can state the algorithm. The transition matrix of the modified Swendsen- 
Wang dynamics (for Qp) on the planar graph G is defined by 

M := M° q = T G>p>q DT GD!r>q T GDjP * >q D*T G;P;q , (9) 

where Tg^, T Gpq and p* are defined in ([2]), Q and ([8]). This dynamics is reversible 
with respect to ir^ q . Although this seems to be complicated, it has an easy inter- 
pretation. The dynamics defined by the transition matrix M performs the following 
steps: 

1) Given a Potts configuration at on G, generate a random cluster state A C Eg 
with respect to TG )Pt q(<r, ■). 

2) Make one step of the Swendsen-Wang dynamics Pp* D q starting at 4o C Eg d 
to get a random cluster state B D C E Gd - 

3) Generate <j t +i with respect to T G (B, •). 

Remark 8. A similar, and even simpler, algorithm was considered by Edwards and 
Sokal [Sokj around 1988 for the two-dimensional torus, which is "self-dual" in some 
sense. In this case they did not have the exact (planar) duality as given in (j7j), but 
one can use a Metropolis acceptance/rejection rule to go to a "dual configuration". 
This algorithm could be advantageous near the critical temperature /3 c {q), but away 
from it the Metropolis step is rejected with probability close to 1. 

The next result shows that the modified Swendsen-Wang has a larger spectral 
gap than the original Swendsen-Wang on G and Gd, respectively. 

Proposition 9. Let P9 be the Swendsen-Wang dynamics on a planar graph G, 
which is reversible with respect to 7r9 , and M as in (Q. Then 

A(M) > max{A(P /3 G 9 ), A(P^)}. 
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Proof. Let tc := TT^ q and n* := as well as \i := \x^ q and // := /i^. Addition- 
ally, let T x := T GiPi9 , T 2 := and 5 := D Tq d p * q Tq d >p * , q D* . Then, since M is 
reversible, for all n e M (recall that S n (a, r) = 7t(t)) 



I /f 1 1 7T 



|(2iST 2 



\{T x BT 2 -S,) n+ \ 

= WT^BTzT^BTz-S^ 



< \\Ti\ 



\B\\t \\T, 



1 (tt,m) 



\{T X BT 2 



,n— 1 



< ims^-s; 



in-l 



l n— 1 



where n := 1 1 Ti | j x 3 (^)-^-x, 2 (w) • ^ ne ^ as ^ inequality holds because the operator 

norms are all less than one (since ttT x = fi, fiB = fi and fiT 2 = n). Taking the n-th 
root and let n tends to infinity we obtain 



\M-SJ 



lim IKSTaTx-^ 



il/n 



< IIS 



< H-B — Sfj,\\ ||T 2 Ti — Sp] 



\B - S„ 



pG _ C 

p,q a* 



with P p G from ©. Since 5 M = D ^» £>* we get 



IM-SUI < 



^p,q M 



< min{||P^-^ 
This yields the result, because A(M) = 1 — \\M — SVII^ 



< II p g d _ q 



| pG _ C 



II P G - 9 



□ 



4.3 Examples 

Now we present lower bounds on the spectral gap of the modified Swendsen-Wang 
dynamics for two examples. The first one is the Potts model on trees, which is mainly 
a toy example, because one knows that the original Swendsen-Wang is rapidly mixing 
at all temperatures |CF98j . We present this example to show that the modified 
Swendsen-Wang can improve the original one at all temperatures. 

The second example is the Potts model on the square lattice, where we prove 
rapid mixing for all non-critical temperatures. 
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4.3.1 Potts model on trees 

Let the graph G be a tree. By definition, G has only one face: the outer one. 
Hence, the dual graph Gd has only one vertex with \Eq\ loops. Therefore the 
Swendsen-Wang dynamics P on Gd is simply independent bond percolation. It 
follows A(P) = 1 and by Proposition we get 

Corollary 10. Let G be a tree and M be the transition matrix of the modified 
Swendsen-Wang dynamics as in Q. Then we get for all > and q G N that 

A(M) = 1. 

4.3.2 The Potts model on the square lattice 

Now we consider the Potts model on the two-dimensional square lattice of side length 
L. This is the graph G = G L = (V, E) with V = {1, . . . , L} 2 and E = {{u, v} G V 2 : 
\u — v\ = 1}. For the definition of the dual graph we use the following notation of 
the union of graphs. Let G = (V, E), W C V and v* V be an additional auxiliary 
vertex. Then we denote by G' = GU w v * the graph G' = (V, E') with V' = VU{v*} 
and E' = E U {{v*,u} : u G W}. Furthermore we denote the "boundary" of Gl by 
h = {(vi,v 2 ) eZ 2 : vi G {1,L} or v 2 G {1, L}} C Vg l . 

By the construction of the dual graph it is easy to obtain (see Figure [TJ that the 
dual graph of the square lattice G = Gl is given by 

G D = Gl-i U Sl _ ± v*. 

We cannot apply Theorem [1] directly for this graph, because it has no bounded 
degree, i.e. deg Go (t>*) = \Sl-i\ = ^{L — 1). Hence, we have to prove the following 
modification of Theorem [TJ 

Theorem [TJ. Let P be the transition matrix of the Swendsen-Wang dynamics, which 
is reversible with respect to TTp q . Furthermore let v G Vq be any vertex, k G [q], 

A k v := {o-Etlp: cr(v) = k} 

and 

„ , , 1 A , ir„( g |A*) V' 

p »<"<t) = —A' + m^) (10) 
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for (X, r G Aj with \a — r\ = 1. 
TZien 



where 



with 



A(P) > cswA(Pl £ ), 
c S w = csw(G,/3,g) := ^O^) 4A 



A := max deg G (w). 

uevb\{«} 



Proof. The proof of this theorem is very similar to the proof of Theorem [U First 
we define the following "flip" transition matrices. Let a G flp and re Aj, then 



Fx(a,r) := fi 



r = a 



a(v) + fc) 



and 



1 9-1 

F 2 (a,r) := F*^) = -^1(t = * + 1) 



1=0 



where (a + Z)(tt) := (er(tt) + I — 1 mod g) + 1 for a G Qp and Z G Z. Let 7r := 7Tg g 
and 7r := ^^(-[A*). It is easy to check that nFi = tt and ttF 2 = n. Following the 
same ideas as in Section |3] with 

Q = F 1 P Akv F 2 PF 1 P Akv F 2l 

which is reversible with respect to 7r = it9 a , we get 



Q S T 



< P 



A 1 ; 



Sir 



P 



AS 



Similarly as Lemma El A(<2) > A(P? fc ). It remains to prove that A(P) > cswA(Q)- 
For this we define for a G Qp, a k := cr — cr(u) + ieAj and a = r {3Z G [g] : 
t = a 1 }. By the construction of the Swendsen-Wang dynamics we have for k, I G [q] 
that 

P(a k ,r l ) = P(a,T) Va,rGfi P . 
Hence we get for a, r G flp with 



c := max — 

CTi,o- 2 ,Ti,T2eAj P(0"2) r 2j 

f71~(T2, T1~T2 
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that 



Q(a,r) = P a £ (^,^i) F 2 {a lt (T 2 ) P(a 2 ,r 2 ) P^r^n) F 2 (n,r) 

<ri,neAj[, 

( T 2 5 T l) 

CTl~CT fe T1=T , 



T2IT2 ~T1 



n = r T2:r 2 fe ~ Tl 



But ri is unique, because t\ = r and n G A„. Therefore, 
Q(ct,t) < cP(a,r) - ]T P^C^n) 



9 

T 2 :T| : ~r fe 



9 

'/ 

< qcP(a,r) 



i=l r 2 eA[:T*v 



With C3 from Lemma [7] we get for o~i, (72, Ti, T2 G with <7i ~ o 2 and n ~ T2 (since 
<J\(y) = cr 2 (v) and T\{v) = t 2 (v)) that 

P(0"l,n) < A P(0"2,Ti) < 2A P(0"2,T 2 ) 2 A 

P(^2,r 2 ) - ° 3 PKr 2 ) - ° 3 P( ( T2,r 2 ) ° 3 ' 
By the same ideas as in the proof of Theorem [I] we conclude 

A(P) > ^A(Q) > ^c^XiQ) > ^(^)~ 4A A(Q). 



This completes the proof. 



□ 



Now we can apply Theorem [TJ for Gp = Gl-i U ( 5 i _ 1 t>* with t> = v* , because 
then A = 4. For the result on the spectral gap of the Markov chain with transition 
matrix M we need a lower bound on the spectral gap of P A k on Go- But this is 

v* 
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simply the heat-bath dynamics on Gl-i with constant (i.e. k) boundary condition. 
From the well-known results on the spectral gap of the heat-bath dynamics on the 
square lattice at high temperatures (see |MQ94] or |MOS94j ) we know that the O(N) 
bound on the inverse-gap also holds for arbitrary boundary conditions. Hence, we 
get the following. 

Corollary 11. Let Gl be the square lattice of side-length L, N = L 2 and M = 

as in OH]) • Then there exist constants eg > and m > such that 

• A(M)" 1 < Cp N for/3? (3 c (q) = ln(l + Jq) 

• A(M)- 1 < cN m forq = 2 and (3 = (5 C {2), 
where c = c S w(G L ,(3 c , 2). 
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